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Abstract. Given a connected (multi)graph G, consisting of V vertices and / 
lines, we consider a class of multidimensional sums of the general form 

Saim, 712, ... ,ni;{Nv}) 



rii— — oono — — oo n/- — — oo 



where the variables qi {i = 1, . . . , I) are real and positive and the variables A'^„ 
(v = 1, . . . ,V) are integer-valued. (5g(ii, n2, . . . , nj; {Ny}) is a function valued 
in {0, 1} which imposes a series of linear constraints among the summation 
variables n^, determined by the topology of the graph G. 

We prove that these sums, which we call Matsubara sums, can be explicitly 
evaluated by applying a G-dependent linear operator Oq to the evaluation 
of the integral obtained from Sg by replacing the discrete variables by 
continuous real variables Xi and replacing the sums by integrals. 



1. Introduction 

Infinite series are ubiquitous in mathematics. In particular, both elementary 
and special functions are either defined in terms of an infinite series or have series 
representations of one sort or another. For example, the Hurwitz zeta function is 
defined by the series 

oo 

n=0 ^ ^' 

for z € CjRez > 1 and q Cz C, q =/= 0, —1, — 2, • • • . Riemann zeta function is a 
special case of Hurwitz zeta function, C(z) = C(z, 1). 

Conversely, confronted with an infinite series, it is always a legitimate pursuit to 
try to evaluate it in terms of known elementary or special functions. Sometimes, 
as an intermediate step, the evaluation of infinite series can be reduced to the 
evaluation of an integral, which may or may not have a closed form evaluation. 
Consider, for instance, Plana's summation formula [Tj, 

(1.2) E fin) = \m + r f{x)dx + ^ H ^-^=^^^dt, 

n=0 ^ 

valid under certain restrictive growth conditions for the function / in the complex 
domain. 



When applied to the series ( 1.1 1, Plana's formula leads to Hermite's representation |! 
. . . . 1 , 1 1 . f°° sin (ztan-i(Vg)) 
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for Reg > 0. This representation actually provides a meromorphic extension of 
C{z,q) to the whole complex z-plane. 

A formula similar in spirit to Plana's summation formula can be obtained by use 
of contour integration and the residue theorem, 



(1.4) 



f^"^^ ^^if 27rnb(z)/(-iz)rfz, 



provided the function ]{z) does not have poles on the imaginary axis. The contour 
denoted above by |t runs parallel to the imaginary axis, upwards from the right 
and downwards from the left, encircling counterclockwise all the poles (located at 
z = m, n S Z) of the kernel nb(2:), defined as 

1 1 



(1.5) 



nb(^) 



(cothTTz — 1) . 



g27rz _ I 2 

If the function f(—iz) is such that its integral along a circular contour at infinity 
vanishes, then we can evaluate the contour integral in (1.4 1 by splitting the original 
contour If into two closed clockwise contours, one on each side of the imaginary axis 
and enclosing all the poles of f{—iz). As a simple example, consider the evaluation 
of the sum 



(1.6) 



E 



1 



where g is a positive real variable. In this case we have the evaluation 

1 



^ (f 27rnb(z) 1^ — ^dz: 



2tt Res nb(2:) — 

z—±q 

= ^ [^b(g) - nb(-g)] 



The identity 

(1.7) nb(z)+nb(-z) + l =0 

can be used to obtain the known result, 

1 TT coth TTq 



(1.8) 



E 



In this paper we consider the evaluation of a class of sums, 971, that generalizes 
the simple sum (1.6). Each of the sums in 971 is defined in terms of a particular 
kind of connected graph, in a way that we make explicit in the next section. 

In principle, the sums in 971, called Matsubara sums, can be evaluated by the re- 
peated direct application of the contour integration formula ( 1.4 1, on a case by case 
basis. However, using an algebraic identity, M. Gaudin [H] has been able to obtain 
a closed form evaluation of any Matsubara sum as a sum of terms corresponding 
to the trees of the corresponding graph. This is reviewed in section [4] 

Starting from Gaudin's result, we prove that any Matsubara sum can be alterna- 
tively evaluated by applying a linear operator to the evaluation of an integral asso- 
ciated with the sum. Although this integral can also be computed using Gaudin's 
method, it is usually the case that the direct computation of the integral can be 
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done in a straightforward manner by other means, including symbohc manipula- 
tion programs such as Mathematica or Maple. In this case the evaluation of the 
corresponding sum will be notably simplified. 

The general results presented in this paper are a by-product of the study of 
general properties of Feynman graphs in the so-called Euclidean or imaginary-time 
formalism of finite-temperature quantum field theory. In this formalism the evalu- 
ation of each Feynman graph requires the computation of a sum of the Matsubara 
type. The existence of the linear operator referred to above was first conjectured 
from the analysis of two general classes of Feynman graphs '6 and then established 
in full generality jJTj . The aim of this paper is to present the relevant mathematical 
results for a readership of non-(particle) physicists. Accordingly, all reference to 
physical quantities has been removed. 

Before defining our class 071 in full generality and presenting the general theorems, 
we will illustrate our main results for the simplest non-trivial sum in DJl, which we 
will call Sg2- 



(1.9) 



OO OO e oo ^ 



^ ^ (,,2+ 2)(„2_^ 2) ^ (n2 + g2)((Ar_„)2+ 2)- 

ni— — oon2 — — oo ^ ^ ^ ' n— — oo ^ ' 

Here e Z and gi > 0, 92 > are real variables. The Kronecker delta symbol is 
defined for n, m G Z by 



(1.10) (5„,„, = 



1 if n = m, 
otherwise. 



A standard evaluation, by the method of residues (|1.4|) for example, yields 
" 1 27r 



(1-11) E 



1 +nb(gi) +nb(g2) 
Ni + qi + q2 



nb(gi) - nb(g2) _ ^hjqi) - nb(g2) _ 1 +nb(gi) + nb(g2) 
Ni - qi+ q2 Ni + qi - q2 Ni - qi- q2 



where n\^{q) is the kernel defined in (1.5 1 



We note that, although written using the imaginary unit i — T, the result 
is clearly real. In the form just written, it becomes easier to recognize the general 
structure of the results we shall present below. 

Consider the integral obtained by replacing, in the LHS of ( |1.11[ ), the sum over 
the discrete variable n by an integral over a continuous variable x: 



Ig2 '■= 
(1.12) 



dx TT{qi + (72) 



(x2 -t- g2) ((TV - x)2 + ql) q,q2 {N^ + {q, + 92)^) 



27r 



2gi2g2 



1 



Ni + qi + q2 Ni - qi- q2 



It is clear from these explicit evaluations, that for this simple case there is a 
simple formal relation between the sum and the integral /gj- To formulate 
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this relation wc introduce reflection operators Ri = R{qi), acting on the space of 
functions of the variables qj, 

(1.13) R{q,)f{qi,...,qi,...) = f{qi,...,~q„...). 

For instance, we have 



and 



Riqi) 



R{q2) 



27r 



2gi2(72 Ni- qi~ q2 



2tt 



2tt 



1 



2qi2q2 iVi + 91 - 92 ' 



2(71 2g2 Ni + qi + q2 



2-K 



1 



2qx2q2 Ni + qi- q2 



In terms of the reflection operator (1.131 we have: 
(1.14) 

Sg, (N, gi, (72) = [1 + nb(gi) (1 - R{qi)) + n^{q2) (l - R(.<l2))] Ig,{N, 9i, 92). 
For instance, the term 

27r nb(gi) 
2qi2q2 Ni + qi - q2 



in (1.11) is generated as 

27r nb(gi) 



2qi2q2 Ni + qi - q2 



= -nb(gi)^(gi) 



2n 



1 



2qi2q2 Ni- qi- q2 



It is straightforward to show that the operator (l — i?i)(l — i?2) annihilates the 
integral evaluated in (1.12), that is: 



(1.15) 



(l-i?l)(l-i?2)/G. =0, 



where we have abbreviated Ri :— R{qi). Therefore, the operator in (1.14) that 
generates the sum from the corresponding integral Ig2 can be written in the 
multiplicative form, 

2 

(1.16) Og, 

where Ubi := nb(gi). 

In this paper we shall prove that all the sums in the class to be defined in 
the next section, satisfy a property similar to (1.14), with an operator of the type 

dTTel). 



Note 1.1. It is important to notice that whereas the evaluation (1.12) is still valid 
if the variable N is extended to the real or complex domains, the same does not 
happen in the case of ( |1.11| . In the latter case one finds, for iV G C, 



(1-17) E 



2tt 



{n^ + ql){{N-nY + ql) 2q,2q2 



l+nb(gi)+nb(g2 + A^^) 
Ni + qi+ 92 



nb(gi) - nb(g2 + A^«) n],{qi) - ii\,{q2 - N i) I + iib{qi) + ii\,{q2 - Ni) 



Ni - qi+ q2 



Ni + qi - q2 



Ni- qi- 92 



Therefore it is clear that the simple relationship (1.14) only holds in the case N £ Z. 
The Matsubara sums to be defined next will in general be functions of several integer 
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variables . We will not consider in this paper the extension of our results to the 
case complex. 



2. MATSUBARA SUMS 

In this section we will introduce the class of sums 9Jl, whose elements we shall call 
Matsubara Sums. These sums were considered for the first time by T. Matsubara 
[3] in his work on the statistical mechanics of quantum fields, where they appear 
in connection to the evaluation of so-called Feynman diagrams. In order to give 
a general definition of the class 371 we need to make use of some graph-theoretical 
terminology. 

Consider a connected graph formed by a set of V points (also called vertices) 
and / edges (also called lines or arcs). We will demand that each line joins two 
different vertices (that is, we exclude loops, i.e. lines that join a vertex to itself) 
and that each vertex be the endpoint of at least two different lines. Any pair of 
vertices can be joined by more than one line. In graph-theoretic language [1][5], we 
are considering a (V, /)-multigraph such that the degree of each vertex is at least 
2. 

We will restrict our attention to graphs of the type described, which we shall 
call Matsubara graphs. 

Let G be a Matsubara graph with V vertices and / lines. We choose for each line 
a definite orientation and assign to this oriented line a positive real number qi and 
an integer- valued summation variable n;. We assign to each vertex v an integer Ny 
and the algebraic sum Ty := SiUi, where 

-1 if the line i is oriented into vertex v 
(2.1) ■Si' = ^ ^1 if the line i is oriented away from vertex v 

if the line i is not incident on vertex v. 

The Matsubara sum Sq associated to the graph G is defined as 

6G{ni,n2,...,nj;{Ny}) 



(2.2) So:= E E 



m — — OO n2— — oo n/— — oo 



(«i + Qt) {nl + ql) ■■■{n] + qj) ' 



Here 5G(ni, n2, . . . , tt./; {A^i,}) is a function valued in {0,1} whose function is to 
impose a series of linear constraints among the summation variables rij. It is given 
explicitly by 



(2.3) 5G{ni,n2,...,ni-{Ny})^\[5T^ 



where 5ra,n is Kronecker delta defined in (1.10 1 above. Whenever there is no possi 



bility of confusion we will use the shorthand notation dcin; N) to denote the object 
5G{ni,n2,...,ni;{Ny}). 

Lemma 2.1. Sg = 0, unless the integers Ny satisfy the relation 

V 

(2.4) Y,Nv^O. 

v=l 
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Figure 1. The graph G2, with labeled vertices and lines (left) and 
its Matsubara sum variables (right). 




Figure 2. The graph G3 with its Matsubara sum variables. 



Proof. The V equations Ty = Ny have to be satisfied simultaneously for the sum 
Sg not to vanish. Now, each summation variable Ui (associated to line i) appears 
in only two of these equations (those corresponding to the vertices on which line i 
is incident), once in the form +ni and once in the form — n^. Therefore, T!„ = 
and 

V V 

a 

Hence, in order for the Matsubara sums to be considered not to vanish identically. 



we shall always assume that the condition (2.4) holds 



We now give a few examples of Matsubara sums. 

Example 2.2. The simplest non-trivial example of a Matsubara sum corresponds 
to the (2, 2)-graph G2 containing 2 vertices joined together by two lines, with iVi + 
N2 = 0, shown in figure [T] Setting A^i = TV, the Matsubara sum for G2 is simply 



the sum Sq^ defined in (1.9) in the Introduction 



Example 2.3. A slightly more complicated example of a Matsubara sum is the 
one associated to the (2, 3)-graph G3, consisting of two vertices joined now by three 
lines, represented in figure[2] Again we require N2 — —Ni = —N, in which case the 
constraints at the vertices reduce to the single equation ni + 712 + = N . Thus, 

(2.5) SGAN,qi,q2,q3) = J2T. ^ 



{nl + ql) (ni + ql) {{N - m - n^)^ + g|) ' 
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Figure 3. The graph Gi, with labeled vertices and lines. The 
orientation of the lines, not shown, is described in the text. 



where it is from now on understood that the summation variables rii run from —oo 
to +00. 

As it was the case for the sum 1 it turns out that Sq^ can also be generated 
in a simple way from the corresponding integral Ig-^ , obtained replacing the double 
sum by a double integral, 

lG3{N,qi,q2,q3) ^ / dxi / dx2-j-^ 



(2.6) 



2gi 2^2 2(73 



■■{xl + qi){xl + ql){{N-x,-X2Y + ql) 
1 1 



_Ni + qi+ q2 + qz iVi - gi - 52 - 93 . 
In this case, after a lengthy evaluation, we find, 

(2.7) S'G3(iV, (71,^2, 93) = OG^{qi,q2,q3)lG3{N,qi,q2,q3), 
where O'q^ is the operator, 

= 1 + nb(<Zi) (1 - Riqi)) + MQ2){1 - R{q2)) + nb(g3) (l - R{.qi}) 
+ nb(<Zi) nb(q2) (l - R{qi)) (l - ^2)) 
+ nb(<Zi) nb(<Z3) (1 - R{qi)) (l - Riq^)) 

(2.8) + nb(q2) nb(q3) (l - ^2)) (l - Riq^)) ■ 



As in the case of G2, it is direct to show that the operator (l— (l— i?2) {^^R-i) 
annihilates : 

(2.9) (l-i?i)(l-i?2)(l-i?3)/G3 =0. 

Therefore, the operator that generates the sum Sg^ from the corresponding integral 
Ig3 can be written in the multiplicative form, 

3 

(2.10) ^G,=\{[l+n^^{l-R^ 



We will prove in the following sections that the types of relationship just de- 
scribed are generic for Matsubara sums. 
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Example 2.4. As a final, more intricate example we consider the Matsubara sum 
associated to the graph G4, represented in figure [s] G4 is a (4,5)-graph. Its 
Matsubara sum will be a function of 3 integer variables, Na, Ni,, Nc, and 5 real 
positive variables, qi, . . . ,q^: 

/r) 1 -. \ o ^ni+n2,Na^n3 — ni—ni,,Ni,^—n3 — ni,Nc 

^' ' „^^„^ K + g?)(ni + <7i)(ni + qDinl + + gf) ' 

The orientation of the lines has been chosen such that line 5 flows from top to 
bottom and all the rest flow from right to left. Solving for the constraints imposed 
by the Kronecker deltas, we flnd that Sd is equivalent to the double sum 

00 ^ ^ 



(2-12) E E { 

n— — 00 m— — 00 ^ 



{n^ + ql){m^ + qi) 



{{Na - n)2 + q^){{Nb + n + mf + ql){{Nb + N, + n + m)2 + qj) 
We shall present the explicit evaluation of this sum in section [6) 

3. The main theorems 

Here we state the main theorems concerning the explicit evaluation of Matsubara 
sums. These results were flrst conjectured [6J and then proved in the context of 
thermal quantum field theory. Our proofs of theorems 3.1 and 3.4 and of lemma 
|3.3| rely on the explicit evaluation of both the Matsubara sum and the Matsubara 
integral associated to a general graph G, which will be given in the next section. 

Theorem 3.1. Let G be a Matsubara graph. Let {qi\,i = 1,...,/ be the set of 

positive real values associated to the lines of G and {-/V^,}, v — 1, . . . ,V be the set of 
integer values associated to the vertices of G. We assume the condition Ny — 
to hold. Then the Matsubara sum ofG, S'cdiVi,}, {(7.;}), can be explicitly evaluated 
as 

(3.1) Sg{{N,}, {gj) = dG{{qi})lG{{N.}. {ft}), 
where 

I 

(3.2) dG{{qi\)^X{[l + n^,,{l-Ri 

i=l 

and Ig{{^v}t {li}) is the Matsubara integral of G. 

Upon expansion of the product, the operator Og can be seen to contain one or 
more terms that individually annihilate the Matsubara integral Ig, so that actually 
the operator Og can be made "smaller" . To state this result we need the following 
graph-theoretical definition [5]: 

Definition 3.2. A cutset of the (connected) graph G is a set of lines whose removal 
from the graph results in a disconnected graph. 

Lemma 3.3. Let C be a cutset of the graph G. Then the operator 

(3.3) = [] (1 _ 4) 

iec 

annihilates the Matsubara integral Ig associated to the graph G. 
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Theorem 3.4. Let the conditions of theorem \3.1\ hold and let L := I — V + 1 be 

the number of independent cycles of the graph G. Then the relationship between 
the Matsubara sum 5g({A'„}, {g;}) and integral lQ({Ny}, {qi}) can alternatively be 
written as 

(3.4) Sg{{N,}, {q,}) = ^'a{{q^})IG{{N,}, {q^}), 

where 

I 

^'a{{q^}) (1 - -^0 + Y! Hb,, (l - R^,) {l - R^,) 



(3.5) 



{iu...,iL) i = l 



Here the indices ii,i2, . . . run from 1 to I (the number of lines of the graph G) 
and the symbol (ii, . . . , ik) stands for an unordered k-tuple with no repeated indices, 
representing a particular set of lines. The prime on the summation symbols imply 
that the tuples that are cutsets of the graph G are to be excluded from the sums. 

Note that the operator O^dgi}) contains products of at most L kernel factors 
nb((7i), since for a graph with L independent cycles the maximum number of lines 
that can be removed without disconnecting the graph is precisely L. 

4. Explicit evaluations of Sq and Iq 
We notice that the main building block of a Matsubara sum can be expressed as 



1 1 



g2 2q 



1 1 



q in — q 



(4.1) = ^h- R{q))—^. 

^ ' 2q^ ^^"q~in 

The representation above would allow us trade the original quadratic denominators 



in the Matsubara sum (2.2 1 for linear denominators (in the summation variables 



rii), and express Sq in the form 

(4.2) Sa^\{—{l-R:^ ^ Mm, n^, . . . , n,; {iVj) 



k=i ^ ^ nf^,n, ~ «"i)(92 - in2) ■ ■ ■ {qi - ini) ' 

were it not for the fact the new sum does not converge in general. However, as we 



shall see, it is possible to regulate the sum in (4.2 1 in such a way that it is well 
defined. 

The main results that will allow us to obtain a complete evaluation of the Mat- 
subara sum Sg were obtained by M. Gaudin [S] long ago, and will be reviewed in 
this section, adapted to the context of this paper. Gaudin showed that the sum- 



mand in (4.2 1 admits a decomposition into partial fractions, which allows us to 
systematically eliminate the constrains imposed by the delta function ^G(n; N) and 
perform the sum explicitly. 



^The number of independent cycles of a graph G is called the cycle rank or cyclomatic number 
in graph theory and is given by m(G) = 7— V + lifGis connected [5]. 
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ni 




Figure 4. The trees of the graph G-^, shown as dark Unas. 

The generalized Kronecker delta Sciji] N) enforces V —1 independent linear rela- 
tions satisfied by the summation variables Ui, also involving the vertex parameters 
N^, which we shall write as 

(4.3) i?„(7V,7i) = 0, for w = 1,...,!/- 1. 

This system of linear equations allows us to solve for F — 1 of the / summation 
variables in terms of a set of L = / — V + 1 independent ones. In general, there will 
be several distinct ways of choosing this set of independent summation variables. As 
shown by GaudinfS^, there is a one-to-one correspondence between the collection 
of all possible sets of independent summation variables and the set of all trees 
associated to the given (connected) graph G. 

A tresis a set of lines of G joining all vertices and making a connected graph 
with no cycles. Every tree T will contain V — 1 lines and its complement T (the 
set of lines of G which do not belong to T) will have L lines. The summation 
variables corresponding to the L lines in T, denoted by n/, will constitute a set of 



independent summation variables in terms of which the system (4.3 1 can be solved. 
The summation variables associated with the lines of the tree, rij, with j G T, 
will be linear combinations of the independent summation variables and the vertex 
parameters, 

(4.4) n,^nJ{N,ni), jeT,leT. 

As a simple example, in figure |4] we show the three possible trees for the (2, 3)- 
graph G^. In this case, each tree T is composed by a single line (heavy line), whose 
summation variable can be expressed, after solving for the constraint at one of the 
vertices imposed by the delta function, in terms of the two independent summation 
variables associated with the (thin) lines that do not belong to the tree (these are 
the lines in T) and the vertex parameter N. For instance, for the first tree we have 
ni = iV — 712 — ^3, etc. 

Similarly, for the (4, 5)-graph G4 of figure [3] there exist 8 trees, shown in figure 
[5] Th e system of equations determined by the Kronecker deltas in the definition 



(2.111 of can be solved in terms of 8 different sets of independent variables. 



one for each of the 8 trees. For instance, the solution that led to the double sum 



(2.12 1, n2 = Na — ni, 713 = Nt + ni + nci and 714 = —Nb — Nc — ni — n^, corresponds 
to the last tree in figure |5] (the one in the bottom-right corner) . 



'^Moro accurately, a spanning tree, in the graph-theoretical terminolo] 



■gy- 
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Figure 5. The 8 trees of the graph G4, represented by the dark hnes. 



in (4.21 



Gaudin's main insight is the following identity for the rational function appearing 



(4.5) t[^^SG{n;N)=J2l[ 



k=l 



Qk - ink 



(4.6) 



For instance, for the case of the Matsubara sum this identity takes the form 
S {N — ni — — "-3) S {N — ni — n2 " ^3) 



((ji - mi) (92 - in2) (93 - ins) {qi + 92 + 93 - iN) 
1 1 



1 



{q2 - in2) (qs - ins) (gi - «"-i) (93 - «"3) (gi - «"-i) (92 - ^?^2) 



We note that the identity (4.5 1 holds only in the presence of the constraint Sain; N) 
on the variables n.^. For instance, in the example above. 



1 



1 



1 



{q2 - in2) {qs - in^) (gi - mi) (53 - mg) {qi - mi) (92 - ^?^2) 

gi + 92 + 93 - »("-i + "2 + n^) 
{qi - mi)(g2 - in2){qs. - in^) ' 



so that the identity (4.6 1 holds only if the constraint ni + n2 + n^ = N is imposed. 



Once we have applied Gaudin's identity (4.51 to the Matsubara sum (4.2 1, we 



proceed to perform the sum over the variables n^. Since, for each tree, the sum- 
mation variables nj^j G T, appear now only in the constraint dG(n;N), we have 
formally 



(4.7) 

Unfortunately, the sum 
(4.8) 



J- J. rij — Ttli * ^ J- A ri, 771. J. A { J rii 1 



ni,...Mi ifzr 



qi - tni ^^^qi-mi 
m-.ieTieT leT 



00 ^ 



q — m 



diverges, and the result above does not make sense. 

However, as shown in [5], it is possible to regulate the sum in such a way that all 
intermediate steps are mathematically sound. The basic idea is to associate to each 
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line i in the graph G a real parameter r^, which at the end will be taken to zero, 
and consider the sum 

gi(niri+n2r2H hnjTj) 

so that 

(4.10) = lim S^. 

The subtle issues related to the interchange in the order of limits and summations 
are discussed at length in |^. Now we have 

1 piniTi 

(4.11) V TT Sain- Ar)e'("i^i+»2^2+-+"^^^) = 17^ ■ 

ni.,...,ni ifzr leT 

where T/ is a linear combination of the r^, whose particular form depends on the tree 
T being considered. It can be shown that Ti is the algebraic sum of the r-variables 
of the lines of the cycle TU {/} (formed by adding the line I to the tree T): will 
be preceded by a plus sign if the line k has the same orientation as the line I, and 
by a minus sign otherwise. 
The sum that we now need is given by 

E inT 
- = 2^enb(eg)e^^ 
q — tn 

n 

where e is the sign of T (e = +1 if T > and e = — 1 if T < 0) and nb(2) is the 



kernel introduced in (1.5). Clearly, only e matters when the regulator T is taken 
to zero. 

The final result for the Matsubara for the graph G is therefore 

(4.13) sa ^ n i^s^ - E (n n 2„,na.,.)^ 

where ei is the sign of the variable T/ associated to each line / e T in the process 
of regulation. 

Consider now the (regulated) Matsubara integral associated to the graph G: 

(4.14) y • • • y dxidx2 ■ ■ ■ dxj 



i{xiTi+X2T2-\ \-XiTi) 



{xl + qf) {xl + qj) ■ ■ ■ {xj + qj) 



Sg{xi,X2, ■■.,xi; {N^}), 



so that 



(4.15) /G=lim/5. 

Now 6g stands for a product of Dirac delta functions, and the integrals over the 
X- variables run from — oo to oo. 
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All the algebraic manipulations given above for the sum Sg hold here unchanged, 
leading us to 

(4.16) [ dxidx2---dxiY\ ^- (5G(a;,iV)e'(^i^i+^^^^+-+^^^^) 

J J ^^qi-ixi 

poo ^ixiTi 

= 11/ —dxi, 

where T; is the same linear combination of the t; as in the case of Sq- 
The integral that we need now is given by 

/oo ixT 
dx = -2TTed(-eq)e^\ 

where e is the sign of T, and d{q) is the Hcaviside step function. Clearly, again 

only £ matters when the regulator T is taken to zero. 

So, 

(4,18) ^ n ^k" - ^ (n ,,-ojW.-.„) n -2-- o^-'"^^ 

5. The relation between Sg and Ig 

In this section we use the explicit evaluations obtained in the previous section to 
prove the theorems stated in section [3] which provide an efficient way of computing 
the Matsubara sum Sg in terms of the Matsubara integral Ig- 

Lemma 5.1. The Matsubara integral can be written as 

(5.1) Ig = Y.^g. 

r 

where 

is the contribution to Ig associated to the tree T. 
Proof. Assuming qi > Q we have 

(1 - Ri) [d{-eiqi)f{qi)] = ^{-eiqi)f{qi) - -disiqi) f {-qi) 
= ~eif{-eiqi). 

Therefore, for each tree T, the action of the operators (l — i?;) with I G T can be 
performed explicitly in (4.181, to yield (5.2 1. Note that ef ~ 1 and that there are 
L = I — V + 1 factors (27r) (one for each line en T) . □ 

Lemma 5.2. For q E R, q ^ 0, the function nb(<z) satisfies the identity 

(5.3) nb((7)--i?(-g) + e(g)nb(k|), 

where e{q) is the sign of q. 



Proof. For q > (5.3 1 is the trivial identity. For q <0, ni,(q) = — 1 — nb(— 9) is an 
immediate consequence of the identity (1.7 1. □ 
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Lemma 5.3. The Matsubara sum Sq can be written as 



(5.4) 



T 



Proof. The identity ( |5.3[ ) implies 

nble/gO = -'&{-eiqi) + e; e(g;) nb(k; |)- 
Then, as we computed already, for qi > and an arbitrary function f{q), 
(1 - Ri) [- d{-eiqi)f{qi)] = eifi-eiqi), 

whereas 

[£/e(gOnb(ki|)/(9i)]=e;nb(|«|) [£(90/(9;) 

= e,e(ffi)nb(|gi|)(l + -R;)/(«)- 



But 



{l + Rl)fi-eiqi) = {l + Ri)fiqi). 



Therefore, for (?/ > 0, 

(1 - ^b(e;9;)/('7/) = l+nbi(l + ^;) f{~£iqi), 



and hence 

n(i-^oin 



27r£; nb(e(9;) 



The final result follows from the basic property. 



(5.5) 



-^{l + R{q))f{q)^{l^R{q)) 



□ 



Now we are finally in a position to prove the propositions of section [3] 
Proof of Theorem\3.1\ We know from lemma [5T3| that 



(5.6) 5G = ^n[l + ^b/(l-i?0]^^ 

But the identity 



r 



(1 - R{q)) i (1 - R{q)) = i (1 + i?(g)) (1 - R{q)) ^ 0, 

implies that all the operators (l — Rj), with j ^ T, annihilate each of the Iq defined 
in (5.2 1. Therefore we can extend the product indices in (5.6) from I G T to i G G: 
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But now observe that the operator acting on Iq is T-independent, so that we can 
transpose it with the sum over trees to get 



ieG T 
= n[l+nb,:(l-4)]/G, 



according to (5.1 1 



□ 



Proof of Lemma \373\ Let C be a cutset of the graph G and let T be an arbitrary 
tree of G. The hues in G cannot all belong to T since then C would not be a cutset 
(recall that T is a connected graph). Therefore, at least one of the lines in G must 
belong to the tree T, say k. But then the operator 



Ac = 11 {I- Rr) 



will contain the factor (l — Rk), which annihilates Iq in (5.2 1. Since this will be 
true for any tree, the result follows. □ 



Proof of Theorem 3.4 The expansion of the product defining the operator Oq gen- 



erates an expression similar to (3.5 1, but where the sums run over all possible tuples 
of lines at each order. But according to lemma [33] if a tuple C is a cutset of G, then 
the corresponding term will contain the operator Ac, defined above, which annihi- 
lates the integral Ic- So the tuples corresponding to cutsets can safely be omitted 
from Oc- Finally, since a tree of G has V — 1 lines, then the maximum number of 
lines that we can remove without disconnecting the graph G is L = I — {V — 1). 
So all tuples with more than L lines will be cutsets and hence the expansion of the 
product defining the operator Oq ends at degree L. □ 

6. Example application: the calculation of 
According to the results presented in this paper, we can compute the sum Sc^ 



defined in (2.11 1 by first computing its associate d in tegral Ic4 and then acting on 
it with the operator Oq^. According to theorem 3.4 this operator is given by 



(6.1) 



l+nbi(l- 


Ri) 


+ nb2 (1 - R2) 






+ nb3 (1 - 




+ nb4 (1 - ^4) + nb5 (1 - 


-R5) 




+ nb2 nb4 


(1- 


i?2)(l-i?4) +nb2nb3(l 


-i?2)(l 


-R3) 


+ nbinb3 


(1- 


i?i)(l-i?3)+nbinb4(l 


-i?i)(l 


-Ri) 


-fnb4nb5 


(1- 


i?4)(l-i?5) +nb3nb5(l 


-i?3)(l 


-R5) 


+ nb2nb5 


(1- 


i?2) (1 - ^5) + nbi nb5 (1 


-i?l)(l 


-R5) 



where, as before, Ub; = nb(<Zi) and Ri = R{qi). 

We note that the operator O'q^ ends at degree 2, since the removal of 3 or more 
lines from G4 disconnects the graph. Moreover, each of the quadratic terms in (6.1 ) 
is associated with one of the trees of G4 shown in figure [5]) . Note that there are no 
quadratic terms in ( |6.1[ ) with index combinations (1,2) and (3,4), since these are 
cutsets of G4 (see figure |3] for the labeHng) . 
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The integral Iq^ associated to G4 is given by (see (2.12 1) 

/oo poo f 

dx dyl 
-00 J —00 L 



1 



1 



am + x + yy + qDUm + N, + x + yr + ql){y^ + ql) , 

This integral can be calculated following Gaudin's approach explained in section |4] 
or by directly performing the x and y integrations, one after the other, as we do 
now. Again, it is convenient to work with linear rather than quadratic denominators 
by expressing 

(6-3) = > 

x'^ + gr^ 2q IX + eq 

so that 

2qi2q22qs2q^2q5 I (^a; + £i<Zi) (i(iV, - x) + £2^2) 

1"°° , £3 £4 £5 

^ ^ {i{Nb + x + y)+ £393) {i{Nb + N^ + x + y)+ £494) («y + £5'?5) 

Using Mathematica 7 to perform first the y-integral and then the x-integral we 
find 



2gi 2^2 2^3 2942^5 



{iNa + qi + q2){l{Na + Nk) + -72 + 9.3 + q5){i{Na + + N^) + ^2 + 94 + ^s) 

1 

" {iNb + qi + q3 + q5){i{Na + Nt) + q2 + qs + q5){i{Na + Nb + N^) + 92 + 94 + 95) 

1 

" (liVa - 91 - q2)iiNb + 91 + 93 + 95)(*(A^fc + A^c) + 9i + 94 + 95) 

1 

' (iTVb + 91 + 93 + q5){i{Nb + N^) + 91 + 94 + 95)(i(A^a + m + N^) + 92 + 94 + 95) 

1 

~ {iNa - 91 - 92)(«iVc - 93 - 94)(i(A^a + Nb + N^) - 92 - 94 - 95) 

1 

(iiV„ + 91 + 92)(«^c - 93 - 94)(«(A^fc + N^) - 91 - 94 - 95) 

1 

^ {iNc - 93 - 94)(i(iVb + A^J - 91 - 94 - 95)(«(A^a + A^h + A^c) - 92 - 94 - 95) 

1 

~ (iAa + 91 + 92)(iAc - 93 - 94)(j(Afa + A^fa) + 92 + 93 + 95) 

1 



{iNa - 91 - q2){iNc - 93 - qi){iNb + 9i + 93 + 95) 
1 

(zA^c - 93 - 94)(«^b + 91 + 93 + q5){i{Na + Nb) + 92 + 93 + 95) 
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{iNa -qi- q2){lNc + qs + g4)(j(^a + Nj,) - 92 - 93 - 95) 

1 

{iNa + qi+ q2){iNc + 93 + q4:){'iNb - qi - qz - qs) 
1 

{iNc + qs + q4:){iNb - gi - 93 - 95)(«(^^a + Nb) - 92 - 93 - 95) 

1 

^ {iNa + 91 + q2){iNc + 93 + 94)(«(A^a + A^b + A^c) + 92 + 94 + 95) 

1 

~ (iiVa - 91 - 92)(iiVc + 93 + 94)(j(A^6 + A^c) + 91 + 94 + 95) 

1 

(iiVc + 93 + q4){i{Nb + N^) + 91 + 94 + 95)(i(A^a + Nb + N^) + 92 + 94 + 95) 

1 

~ {iNa - 91 - 92)(i(iVa + Nb) - 92 - 93 ~ 95)(i(A^a + iVfc + N^) - 92 - 94 - 95) 

1 

(iAffi - 91 - 93 - 95)(*(A^a + A^b) - 92 - 93 - 95)(«(A^a + A^b + iVc) - 92 - 94 " 95) 

1 

{iNa + 91 + q2){iNb - 91 - 93 - qb){i{Nb + N^) - 91 - 94 - 95) 

1 

' (iTVf, - 91 - 93 - 95)(*(A^b + Nr) - 91 - 94 - 95)(«(A^a + A^b + A^c) - 92 - 94 - 95). 
The explicit evaluation of sum Sq^ can be obtained from the application of the 



operator O'q given by (6.1) to the result (6.4) for the integral Iq^- The resulting 



expression would fill several pages of this journal, but it can be easily be generated 
by a symbolic manipulation program such as Mathematica. 
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